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ABSTRACT 

Small-scale dark matter structure within the Milky Way is expected to affect pulsar 
timing. The change in gravitational potential induced by a dark matter halo passing 
near the line of sight to a pulsar would produce a varying delay in the light travel time 
of photons from the pulsar. Individual transits produce an effect that would either 
be too rare or too weak to be detected in 30-year pulsar observations. However, a 
population of dark matter subhalos would be expected to produce a detectable effect 
on the measured properties of pulsars if the subhalos constitute a signihcant fraction 
of the total halo mass. The effect is to increase the dispersion of measured period 
derivatives across the pulsar population. By statistical analysis of the ATNF pulsar 
catalogue, we place an upper limit on this dispersion of logCp < —17.05. We use this 
to place strong upper limits on the number density of ultracompact minihalos within 
the Milky Way. These limits are completely independent of the particle nature of dark 
matter. 
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1 INTRODUCTION 

Dark matter structure is thought to have been seeded by ran¬ 
dom density perturbations in the early universe, collapsing 
well after recombination into the first gravitationally bound 
dark matter halos. These small scale structures are expected 
to have hierarchically merged into the larger structures we 
see today (see e.g. Lemoine et al. 2008 and Mo et al. 2010 
on inflationary cosmology and structure formation, respec¬ 
tively). From CMB observations the spectrum of these pri¬ 
mordial perturbations is expected to be nearly scale-free; the 
power is nearly equal on all spatial scales (Hinshaw et al. 
2013; Planck Collaboration 2014). However, current lim¬ 
its only strongly constrain these fluctuations on very large 
scales (wavenumbers fc < 3 Mpc“^). On much smaller scales, 
many cosmological theories beyond the standard model pre¬ 
dict that the power spectrum may deviate significantly from 
the simple Harrison-Zel’dovich (scale-free) model (Adams 
et al. 1997, 2001; Ashoorioon et al. 2009; Erickcek & Sig- 
urdson 2011). 

If significant additional power is present on smaller 
scales, more fluctuations of very large amplitude {5 > 0.3) 
will be produced than otherwise. Such large-amplitude fluc¬ 
tuations rapidly collapse to form primordial black holes 
(PBHs; Carr & Hawking 1974; Carr 1975). While there are 
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tight constraints on the abundance of these rare objects, it 
was recently proposed that smaller amplitude fluctuations 
(0.3 > (5 > 10“^) may give rise to dense dark matter struc¬ 
tures known as ultracompact minihalos (UCMHs; Berezin¬ 
sky et al. 2003, 2012, 2013; Ricotti & Gould 2009; Scott & 
Sivertsson 2009). Instead of collapsing directly as a black 
hole, these perturbations grow by gravitational accretion 
like any other density perturbation until they collapse: loga¬ 
rithmically during radiation domination, and linearly during 
matter domination. Unlike regular perturbations, the large 
initial value of the overdensity means that they enter the 
non-linear regime of growth (i.e. collapse) far earlier than 
do regular 5 ~ 10“® perturbations responsible for most of 
the large-scale structure seen today. Due to this early time 
of collapse, the infall of dark matter onto UCMHs is es¬ 
sentially radial, and their structure will have a very steep 
density profile as a result (p oc Bertschinger 1985; Ri¬ 

cotti & Gould 2009). Although UCMHs continue to accrete 
both dark matter and baryonic matter after recombination 
up until the current era, this steep profile means that they 
are not expected to be tidally disrupted during the course of 
their evolution (Berezinsky et al. 2006, 2008) - to the extent 
that the probability of survival through to the modern era 
is essentially unity for all UCMHs that we consider in this 
paper. 

The present-day mass of a UCMH or PBH may be di¬ 
rectly mapped to the wavenumber of the primordial fluctua- 
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tion that originally seeded it, as the wavenumber of a mode 
re-entering the horizon at any time after inflation depends 
on the horizon size, which sets the initial mass of the over- 
dense region. Because they require a smaller initial overden¬ 
sity than PBHs, UCMHs are expected to be produced in far 
greater numbers than PBHs. UCMHs are therefore a very 
promising direct link to the conditions of the early universe. 
By constraining their abundance, we can place limits upon 
the conditions that would lead to their formation (e.g. Josan 
& Green 2010; Bringmann et al. 2012; Shandera et al. 2013; 
Berezinsky et al. 2011; Anthonisen et al. 2015). 

Non-detection of PBHs has been used to place very 
weak limits on primordial curvature perturbations, of 
log'PTj < —1.5 for wavenumbers in the range 10“^ ^ k < 
10^® Mpc“^ (Josan et al. 2009; Carr et al. 2010; Alabidi et al. 
2012). Non-detection of UCMHs by gamma-ray searches has 
been used to constrain their present-day number density, im¬ 
proving this limit significantly, but over a far smaller range 
of scales: logV-jz < —6.5 for wavenumbers in the range 
3 < fc < 10^ Mpc“^ (Bringmann et al. 2012). It is important 
to note, however, that the limits from gamma-ray searches 
depend entirely on the assumption that dark matter can 
self-annihilate; indeed, UCMHs have been studied exten¬ 
sively as targets for indirect detection of dark matter due 
to their extremely steep density profiles (Scott & Sivertsson 
2009; Lacki & Beacom 2010; Yang et al. 2011a,b,c, 2013a,b,c; 
Zhang 2011; Zheng et al. 2014). 

Gravitational leasing has long been used as an impor¬ 
tant tool in the detection of dark matter, as gravity is the 
only force that we are certain dark matter interacts via. Non¬ 
observation of characteristic changes in the position or light 
curve of a star due to intervening masses (known respectively 
as astrometric and photometric microlensing) have been 
used or proposed as means to weakly constrain the abun¬ 
dance of both PBHs (Tisserand et al. 2007; Wyrzykowski 
et al. 2011a,b; Griest et al. 2013) and UCMHs (Ricotti & 
Gould 2009; Li et al. 2012; Zackrisson et al. 2013). However, 
a confirmed detection of an individual object would require 
either extreme sensitivity or high abundance, so it is not 
expected that these methods will be able to effectively con¬ 
strain the properties of the primordial power spectrum in 
the near future (unless something like the proposed THEIA 
satellite mission flies). 

As a more sensitive alternative, it has been proposed 
that substructure could be detected by measuring the effect 
of an intervening mass on the timing of a millisecond pul¬ 
sar (Siegel et al. 2007). This ‘time-delay leasing’ uses the 
increased travel time of a light ray that passes through a 
changing gravitational potential, known as the Shapiro ef¬ 
fect. Although it is not possible to directly measure the de¬ 
lay due to a static mass, a dark matter halo that moves 
between an observer and a pulsar would cause the pulse 
frequency to appear to decrease as the lens travels toward 
the line of sight, and increase again as the lens moves away. 
Even though this is a very weak effect, millisecond pulsars 
can provide extremely accurate clocks - in many cases, sig¬ 
nificantly more accurate than atomic clocks - so it may be 
expected that large mass subhalos would indeed produce a 
detectable effect on pulsar timing. 

Here, we use time-delay leasing to derive improved lim¬ 
its on the number density of UCMHs within the Milky Way. 
We first describe the analytical method for calculating the 


time delay produced by a UCMH (Section 2.1). We then 
compute the probability for a dark matter halo to transit 
the line of sight to a pulsar, and to observe the event by the 
delay of the pulsed emission (Section 2.2). We predict the 
population impact of this effect to be a type of Gaussian 
noise present in all pulsar period derivative (P) measure¬ 
ments, and describe a novel method for using this noise to 
constrain the properties of dark matter substructure (Sec¬ 
tion 2.3). We use the resulting limits on P-noise to set lim¬ 
its on the number density of UCMHs within the Milky Way 
(Section 3), then conclude and summarize (Section 4). 


2 METHOD 

2.1 Time-Delay Lensing 

The Shapiro time delay of a light ray may be investigated 
by considering the travel time of the light from its source. 
In the case of an unlensed system, this is simply the proper 
distance divided by the speed of light. However, in the lensed 
case, this will differ by 

Afltt = Atgeo -f Afpot, (1) 

where Atgeo is the geometric component of the delay, i.e. the 
change in the light ray’s path length due to lensing. Atpot 
is the gravitational potential contribution, which may be 
calculated from the Newtonian potential through which the 
light ray passes, as an integral along its path (Betters et al. 
2012) 



c 


Here c is the speed of light, ip is the Newtonian potential, r 
is the radius from the centre of the halo, and C is the path 
of the light beam, parametrized by ds. 

In the case of lensing of sources within the Milky Way 
by small scale dark matter halos, the deflection of the light 
ray will be negligible, and so the difference between a lensed 
system and an unlensed one will be significantly greater in 
the potential term than the geometric. We may therefore 
approximate the path of the light ray as a straight line; we 
refer to this as the zero deflection approximation. The line 
integral may then be solved by setting r = y/s^AAA, so that 


2 

Afitt = Atpot = -? / V’(s) ds, (3) 

■'Dda 

where we define the position of the lens along the path as 
s = 0. Here, b is the impact parameter of the beam, Dd is 
the distance from observer to lens, and Dda is the distance 
from lens to source. 

The Newtonian gravitational potential for any extended 
spherically symmetric mass with radially varying density is 



( 4 ) 


Under the approximation of zero deflection, the light travel 
time from source to observer may be found if the Newtonian 
potential is easily integrable. Here we are specifically seeking 
to investigate dark matter substructure that that has not 
been amenable to detection by standard lensing effects like 
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source magnification or shear, both of which rely on non- 
negligible light deflection. The assumption of zero deflection 
is hence not only mathematically convenient in this case, 
but is implied by the problem itself. 


2.1.1 Ultracompact Minihalos 


UCMHs have been predicted to have density profile (Ricotti 
et al. 2008; Ricotti & Gould 2009) 


p(r, z) = 
where 
l^siz) = 


Ks{z) 
r9/4 ’ 


SfxMhjz) 

16TvRh{z)i 


(5) 

( 6 ) 


Here is the ratio of dark matter to matter in 

the universe, Mh is the mass of the halo at redshift a, and 
Rh is the radius of the halo. 


^ = 0.019 

pc V^ + 1/ V Mq ) 


(7) 


The approximation of purely radial infall breaks down 
at small r, so we take the halo to have a flattened density 
profile within a core of radius Vc. at zero redshift (Bringmann 
et al. 2012), given by 


Tc 

K 


2.9 X 10”'^ 


f 1000 

V 2c + 1 


f M^\ 


( 8 ) 


where Zc is the redshift at which the halo collapsed, 
is the present day mass, and R® is the present day halo 
radius. In the case of annihilating dark matter, the density 
is truncated at a slightly lower value, increasing the radius 
of the core, and decreasing the total mass of the halo by 
less than a percent - having a negligible effect on the time 
delay. Therefore, for the remainder of this paper we consider 
only the density truncation at r^. The particle nature of 
dark matter can also impact the possible masses of UCMHs 
formed in the early Universe, making any limits we draw 
invalid below a certain cutoff mass. This mass can vary from 
10“^ to 1O“^^M0 for typical WIMPs (Bringmann 2009). 

From Eqs. 4 and 5, we find the gravitational potential 
induced by a UCMH to be 


ip{Rh < r) = 


2'kGks I 
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24 \ 
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ip{rc <r < Rh) = 


(9) 


- Mm h (10) 


( 11 ) 


By integrating these potentials along the line of sight 
(Eq. 3) under the assumption of zero deflection, and defin¬ 
ing the lens position along the LOS as s = 0, we find the 
potential component of the light travel time passing from si 
to S 2 within each region of the halo, as shown in Fig. 1. 


Source 



Figure 1. Schematic diagram of undeflected leasing by an ultra¬ 
compact minihalo. The photon path is shown as a dashed line, 
passing from source, through the different regions of the lens, to 
observer, with impact parameter b. Size of core and halo are ex¬ 
aggerated for visibility. 


ii) Within the halo: 


thaloUx ^ S2 ) 


StvGks 12 (s 2 — Si) 
3c3 

I 


16 [ 52 ^( 52 ) - siF’(si)] 

& 1/4 


-I- 3rf 7“* In 


+ sij 


(13) 


where F{s) = 2 F 1 is a Gaussian hypergeo¬ 

metric function. 


iii) Within the core: 


AtvGks , . 

= (Sl, S2) = -^-(S2 - Sl) 
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1/4 
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(3b -t- Sl -fi S 1 S 2 -t- S 2 ) 
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(14) 


i) Outside of the halo: 


text(si, S2) 


StvGks 

3c3 




'Ub^ + s| -b S2\ 
\A^”+^ + Sl / ’ 


( 12 ) 


The total delay may then be found as the sum of each 
section through which the light ray passes. For example, the 
total potential delay of a light ray originating outside of the 
halo, passing through it (without intersecting the core), and 
received by an observer on the exterior of the halo may be 
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found by: 

^pot — ^ext ^-^ds; \J 

+ tnaio [sjRl-b\-^Rl-h^'^ 

+ U^t (-\lRl-h'^,-D^ . (15) 


2.1.2 Navarro-Frenk-White Subhalos 


Similarly, we can also calculate the light travel time for the 
commonly considered Navarro-Frenk-White (NFW) dark 
matter density profile. NFW halos have density profile 


p(r) = 


hcPc 


r/rs (1 -h r/rs) 


5 ’ 


where 

, 200 
dc — 


3 ln(l-I-c) — c/(l-I-c) ’ 


(16) 


(17) 


pc = /SttG is the critical density for closure of the uni¬ 
verse at redshift z = 0, Hq is the present value of the Hubble 
constant, G is the gravitational constant and = r 2 oo/c is 
the scale radius. Here c is the concentration parameter and 
r ’200 is the virial radius, taken as the radius beyond which 
the halo is truncated as p{r > r 2 oo) = 0. 

We adopt the fitting function for the mass- 
concentration relation for NFW halos in the low-redshift 
regime given by Correa et al. (2015), allowing each halo to 
be described solely by its mass, as 


rsc{Mh) = r2oo = 


( SMh 

\ SOOttPc 


1/3 


The Newtonian gravitational potential of such a halo is 


ipir < r 2 oo) = -^TfGScPcrs 


In 


r + Vs 


1 


1 + c 


</>(r 2 oo < r) = - 


A-kGScPcTs 


ln(l -I- c) - 


1 -I- c 


(18) 

s 

(19) 

( 20 ) 


This potential is not analytically integrable, so we integrate 
it numerically in order to compute the potential time delay. 


2.2 Probability of Individual Subhalo Detection 

Although it is not possible to directly measure a light ray’s 
lensed travel time from a source, changes in the leasing sys¬ 
tem will cause a change in the light travel time, which in turn 
will be measured as a change in the residual of the source 
pulsar’s timing. Following the method outlined in section 
2 .1, we may iteratively calculate the time delay as a lens 
moves across the line of sight. We give an example of this 
signal for both NFW and UCMH dark matter profiles, com¬ 
pared to that due to a point mass, in Fig. 2. Due to its 
significantly steeper gravitational potential, a dark matter 
subhalo with a UCMH profile would consistently produce a 
stronger time-delay signal than one following an NFW pro¬ 
file. Due to their steep density profiles, UCMHs are - for 
time delay purposes - essentially pointlike, particularly at 
masses < 1 Mq. This makes them prime candidates for any 



Figure 2. An example of a time-delay signal from halos with both 
UCMH and NFW profiles, as they would be seen in a pulsar’s 
timing residual over a typical pulsar observation time of 30 years, 
selected such that it passes its point of closest approach to the line 
of sight after exactly 15 years. Here we have chosen a halo mass 
Mf^ = 1 X 10® M 0 , impact parameter b = 10 pc, a halo velocity 
perpendicular to the line of sight v± = 200 kms“®, observer- 
pulsar distance Ds = Ikpc, observer-lens distance = 0.5 kpc, 
and for the NFW profile, concentration parameter c = 16. 


gravitational search for dark matter structure (or alterna¬ 
tively, provide strong opportunities to limit their properties 
should none be detected). 

The rate at which detectable transit events occur for 
a given pulsar is entirely dependent upon both the number 
density of halos of a given mass, as well as their velocity 
distribution. This rate may be predicted by construction of 
a simple simulation, in which halos are stochastically dis¬ 
tributed along the line of sight. Here, we approximate the 
local dark matter density to be given by a global NFW pro¬ 
file at r = 8 kpc, with Mvir = 9.4 x c = 18, provid¬ 

ing a local dark matter density of 0.285 GeV.cm“®. We take 
a fraction / of the local density to be confined within halos 
of mass Mh. We assign each of halo a speed of 200 kms“® 
in the Galactic rest frame, and distribute halos’ directions 
of motion isotropically. We likewise assume that both the 
observer and a nearby pulsar (which we set at a distance of 
2 kpc) co-rotate the galactic centre at 220 kms“®, neglecting 
the motion of the Earth around the Sun. Allowing each ob¬ 
ject to continue along its trajectory for 30 years, we record 
the number of halos that pass within a given distance rmax 
of the line of sight to a pulsar, thereby calculating the ap¬ 
proximate rate as a function of rmax for a given fraction / 
and halo mass Mh. 

As each of these transit events occur independently of 
one another, their occurrence can be modelled as a homo¬ 
geneous Poisson process. The probability that at least one 
event will occur in a pulsar’s signal, within a given Tmax and 
observation time r, is 

P>1 transit (t,w) = 1 - (21 ) 
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where A is the event rate determined by simulation. We de¬ 
mand that such an event occur at least once within the 30- 
year observation data of the 315 pulsars present in the ATNF 
catalogue (Manchester et al. 2005) at approximately 2kpc 
from Earth (at >95% CL). These requirements correspond 
to a minimum rate of A > 1.0 x lO”'* yr~^. 

Following the method we describe in Section 2.1 and 
taking the optimal case (in which a halo transits along a 
path perpendicular to the line of sight), we can compute the 
30-year amplitude expected from a transit at radius rmax- 
Additionally, we simplistically assume that if the 30-year 
signal amplitude exceeds some sensitivity threshold, the sig¬ 
nal may then be taken to be detectable within the pulsar 
data. Using Brent’s root-finding algorithm (Brent 1973), we 
vary / for each Mh to find the required event rate, while still 
producing a signal that exceeds the sensitivity threshold. We 
show the substructure properties that would be expected to 
produce at least one detectable event in Fig. 3, for two dif¬ 
ferent assumed timing sensitivities (Ins and 10ns). 

For all NFW masses that we investigated, the halo 
fraction that would produce detectable events remains en¬ 
tirely outside of physically reasonable scenarios (i.e. / > 1). 
UCMHs of mass ~ 10~^Mq would be expected to be seen 
in pulsar timing data if they constitute > 5% of the local 
dark matter. It should be noted, however, that observation 
of nearby pulsars would be expected to only provide limits 
on local dark matter substructure properties. For NFW ha¬ 
los, this is only indicative of the local dumpiness - however, 
UCMHs are immune to tidal disruption and so their number 
density is expected to follow the global NFW profile of the 
Milky Way, and thus the local value of / is expected to be 
the same as the global /mw. 

While for some pulsars the time of arrival of a pulse 
may known with < 50 ns precision (Hobbs 2014), others 
achieve far less (~ l^s). Our assumed sensitivity threshold 
is therefore very optimistic - with present pulsar data, a de¬ 
tection of any dark matter halo would not be expected. In 
the future, if a complete search through timing observations 
of nearby pulsars were to be undertaken at our assumed im¬ 
proved sensitivity, and a timing event found, the discovered 
object could be identihed as an ultracompact minihalo, a 
compact object (baryonic or otherwise), or as a yet unde¬ 
scribed dark matter structure. Even in the most optimistic 
case, in which the entire galactic halo is taken to consist en¬ 
tirely of subhalos of equal mass (/ = 1), we End that if a time 
delay event were to be identified within 30 year data, such 
an object would not be an NFW halo (at > 99.98% CL), un¬ 
less timing accuracy significantly exceeded 1 ns. The positive 
identification of any dark matter substructure, while highly 
difficult, would have significant implications for our under¬ 
standing of the conditions present in the early universe, as 
well as the identity and structure formation of dark matter. 


2.3 Effect on Observed Pulsar Period Derivatives 

Searching through pulsar timing data for individual halo 
transit signals would be an intensive process, and, if a po¬ 
tential detection were found, it would be extremely difficult 
to differentiate from effects like timing noise, glitches and 
polynomial fitting errors in periods and period derivatives. 
A candidate detection could only be confirmed if the signal 
were seen in two or more pulsars near to each other on the 
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Figure 3. The lower limits on the halo mass distribution for 
which > 1 detectable time-delay leasing event would be present 
in the timing signal of the 315 pulsars at distance ~ 2kpc in 
the ATNF catalogue (Manchester et al. 2005), at 95% CL. We 
show limits for both NFW and UCMH halo profiles, and for 30- 
year signal amplitude sensitivities of 1 X 10“® s and 1 X 10“® s. 
We adopted the local dark matter density throughout the 2 kpc 
sphere around the Sun, given by the value at r = 8 kpc in a Milky 
Way NFW profile with c = 18 and = 9.4 X 10*^*^ Mq, again 

corresponding to a local dark matter density of 0.285 GeV/cc 
(Battaglia et al. 2005, 2006). The region in which / > 1 represents 
scenarios that are unphysical, i.e. where the integrated mass of 
all subhalos would exceed the known mass of the Milky Way. 


sky. As we showed in the previous Section, the alignment 
required to produce a transit effect within the 30 year data 
is extremely rare anyway, unless the number density of halos 
is very high. 

Alternatively, instead of looking for individual subhalos, 
searching for net effects on observations of the entire pulsar 
population due to substructure could be a more effective 
way of constraining the dumpiness of the Milky Way’s dark 
matter halo. For every halo that completes a transit across 
the line of sight within 30 years, there are a multitude more 
that are simply travelling toward or away from their point 
of closest approach, and thus contribute to the net signal 
linearly. That is, it is significantly more likely that the total 
observed timing effect will be a sum of linear effects, rather 
than follow the characteristic peaked shape of an individual 
transit shown in Fig. 2. This linear term would ultimately be 
measured as a contribution to the observed period derivative 
of a pulsar Fobs, 

i) r> 1 /^^pot,i(i + At) — tpot,i(t) 

obs ~ Jt^pulsar \ / ^ I I ? 

for N subhalos, calculated over a time interval At. 

The second term of Eq. 22, which we refer to as the 
‘P-noise’, is the sum of a function of random variables. By 
the central limit theorem, it would therefore present as a 
Gaussian with mean of zero in the case of large N, and vari- 
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P (sec) 


Figure 4. The period-period derivative (P-P) distribution of non¬ 
binary pulsars within the Milky Way from the ATNF pulsar cat¬ 
alogue (Manchester et al. 2005), co-plotted with a histogram of 
period derivatives. The upper limit on the standard deviation of 
the Gaussian P-noise due to intervening dark matter is displayed 
as the dashed black line. Binary pulsars have been excluded in or¬ 
der to eliminate potential additional contributions to the period 
derivative. 


ance dependent upon the underlying variables: halo velocity, 
position, mass, and number density. If the standard devia¬ 
tion of the P-noise is comparable to the value of the actual 
period derivative, this additional term would typically dom¬ 
inate over the ‘true’ period derivative of the pulsar. 

By examining the distribution of observed pulsar pe¬ 
riod derivatives from the ATNF catalogue (Manchester et al. 
2005), shown in Fig. 4, it can be seen that the majority 
of non-binary pulsars are observed to have positive period 
derivatives. These appear to follow a log-normal distribu¬ 
tion. It is therefore possible to determine how large the stan¬ 
dard deviation of the Gaussian P-noise may be, such that 
the observed distribution can still be reproduced. For ex¬ 
ample, if ap > 10”^^, the noise would dominate the ‘true’ 
distribution of pulsars, such that they would be observed as 
being strongly normal, rather than log-normal. 

In order to investigate the upper limit on up, we model 
the ‘true’ distribution of pulsar period derivatives as a log¬ 
normal distribution. Given that the luminosity distribu¬ 
tion of isolated pulsars has been shown to be log-normal 
(Faucher-Giguere & Kaspi 2006; Ridley & Lorimer 2010), 
this is a reasonable assumption. This relationship can be un¬ 
derstood by pulsars only losing rotational energy at a rate 
linearly proportional to their luminosity. 

When observed, the true P distribution will appear con¬ 
volved with both the instrumental noise and the time-delay 
noise, 

fobs{P) = exp [N{^, CTtrue)] * A^(0, CTp) * iV(0, uf ) (23) 

where N{fi, is the normal distribution, ap is the induced 
variance due to time delay lensing, ai is the variance due to 
instrumental noise, and a^^^^ are respectively the mean 
and variance of the ‘true’ log-normal distribution of pulsar 


period derivatives. This convolution results in the likelihood 
function for N pulsars of (Hawkins 1991) 


P(Pobs|u'p, fJtrue) — 


)-n/ 

JO 


l-lo 27 r(Ttruel/y^O-p -|- CTj ; 


X exp 


(Iny-g)^ 

2cr?rue 2(0-2 ^2) 


dy. 


(24) 


To determine the allowed values of ap, we fitted the 
observed values of Fobs from the ATNF catalogue using Eq. 
24. Using the publicly-available nested sampling algorithm 
MultiNest v3.9 (Feroz & Hobson 2008; Feroz et al. 2009, 
2013), we scanned over log up, ^ and Utrue with uniform 
priors. We adopted the values of ai given for each individual 
pulsar by the ATNF catalogue. 

We found the posterior mean of the P-noise standard 
deviation, which we obtained by marginalising over the log¬ 
normal parameters ^ and a^me, as logUp = —17.21 - with 
95% credible interval upper limit of —17.05, and lower limit 
of —17.31. However, millisecond pulsars are not so sim¬ 
ply described, and may exhibit a range of complex tim¬ 
ing phenomena contributing to their final measured sig¬ 
nal. It therefore should not be assumed that their spread 
is solely due to intervening dark matter and instrumental 
error. As such, we adopt the upper limit implied by the 95% 
confidence Bayesian credible interval on logUp. Our limit, 
logcTp < —17.05, is shown as the black dashed line in Fig. 4. 
Should the dark matter distribution be such that the P-noise 
is equal to this value, it would be expected that approxi¬ 
mately half of the pulsars with period derivative beneath 
this will be observed to have a negative period derivative. 


3 LIMITS ON ULTRACOMPACT MINIHALO 

NUMBER DENSITY 

We can compare the magnitude of P-noise expected from a 
population of UCMHs with a given mass to the upper limit 
that we found in Section 2.3, using the method outlined in 
Section 2.1 and Eq. 22. To this end, we simulated UCMHs 
distributed stochastically according to a Milky Way-like 
NFW profile, again adopting c = 18, M„ir = 9.4 x IO^^Mq 
from Battaglia et al. (2005, 2006), in order to produce a 
simulated period derivative contribution. 

Galactic dark matter subhalo velocities are currently 
poorly understood, and so we took a simplified model in 
which all subhalos have isotropic velocities of magnitude 
200kms“^, independent of galactocentric radius. We placed 
the observer at a galactocentric radius of 8kpc from the 
centre of the Galactic profile, and simulated observations 
of mock lines of sight to all 1810 non-binary pulsars with 
known period derivatives in the ATNF catalogue. 

As the time delay due to each individual subhalo is addi¬ 
tive and independent, it is expected that the contributions of 
lenses of different masses would likewise be additive. In this 
way, we constrain the fraction of the Milky Way in UCMHs 
of each mass independently. It should therefore be noted 
that models predicting a range of UCMH masses must have 
their predictions integrated over our single-mass limits in 
order to properly assess their validity. 

Varying the fraction of dark matter contained in 
UCMHs of a given mass, and matching the output P-noise 
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to our upper limit from pulsar data, we obtain an upper 
limit on UCMH number density within the Milky Way as a 
function of mass. We show this limit in Fig. 5, comparing to 
those from gamma-ray searches by Bringmann et al. (2012), 
which necessarily assume a specific model for annihilating 
dark matter. Mirroring their procedure, we reduce our limits 
at large masses (> 10^Mq), corresponding to case in which 
the P-noise does exceed the observational upper limit, but 
the large-N condition of the central limit theorem no longer 
holds. By performing a chi-squared test for normality on 
the time delay distribution at low N, we find that that for 
> 12 the distribution remains normal (p-value < 0.05) for 
all halo masses that we consider. As a conservative measure, 
we therefore reduce our limit to > 20 halos within the 
Milky Way, as can be seen in Fig. 5. 

Our final limits on UCMH number density are of simi¬ 
lar magnitude to those found by gamma-ray searches, with 
strongest constraint / < 1.1 x 10“® at a mass of ap¬ 
proximately 1 X lO^M©. In fact, our constraints from pul¬ 
sar timing are significantly stronger throughout the range 
10“^ < Mh/Mq < 10^ than those found previously. Not only 
are these the strongest limits on the abundance of UCMH 
number density to date, they have the additional benefit of 
being the only strong constraint that does not rely on the 
assumption that dark matter undergoes annihilation. 

Baryonic matter along the line of sight may be expected 
to produce the same gravitational effects. While additional 
structure (baryonic or otherwise) would boost this signal 
beyond that provided by UCMHs alone, we have provided 
only upper limits on the total Gaussian noise, and so our 
constraint on UCMH abundance still holds true. In addition 
to gravitational effects, baryonic matter is known to change 
the speed of light propagation - providing an additional (po¬ 
tentially varying) delay. The strength of this effect changes 
on far shorter timescales than that due to dark matter struc¬ 
ture, and so does not change the measured period derivative 
of the pulsar (You et al. 2007). While some pulsars do in¬ 
deed appear to have a constantly increasing or decreasing 
dispersion measure, this would have the same effect on the 
period derivative distribution (at a single wavelength) as the 
effect we investigate - adding to the total Gaussian noise. 


4 CONCLUSIONS 

The difference in photon travel time due to an intervening 
source of gravitational potential has long been known, most 
famously as the Shapiro delay. It was recently proposed that 
movement of dark matter structures along the line of sight 
to a pulsar could potentially produce a measurable change 
in the pulsar’s timing signal. If measurement of these in¬ 
duced shifts were possible, the properties of the intervening 
structures could potentially be investigated, providing in¬ 
sight into the nature of dark matter and the formation of its 
structure. 

UCMHs are one such predicted form of dark matter sub¬ 
structure, expected to be seeded in the very early universe. 
As a consequence of their early formation, they have ex¬ 
tremely steep density profiles, and consequently steep grav¬ 
itational potentials, allowing them to provide a strong time- 
delay signal. UCMHs have been shown to be highly immune 
to disruption by tidal forces, and so are expected to have 


Mh (Me) 



Figure 5. The 95% credible interval upper limits on the fraction 
of dark matter in the Milky Way contained in UCMHs. Limits 
are shown from both pulsar period derivatives and non-detection 
in gamma-rays by Fermi-hAT searches for dark matter annihila¬ 
tion. We map UCMH masses to wavenumber k at horizon entry 
following the method outlined in Bringmann et al. (2012). 


persisted from their collapse shortly after matter-radiation 
equality through to the present day. They therefore provide 
a unique probe of the conditions of the early universe. 

Here we have calculated the probability of detecting an 
individual halo by pulsar timing. We find that detections 
are likely impossible for NFW minihalos, but UCMHs may 
produce a detectable signal, should the number density of 
lenses be appropriately high (/ > 0.01) for UCMHs in the 
mass range 10“^ < M^/Mq < 10^. 

More excitingly, we predict an additive time-delay 
‘noise’ on pulsar period derivatives due to a population of 
dark matter halos. By determining an observational upper 
limit of logCTp < —17.05 on the observed amount of this 
noise, we have placed upper limits on the number densities 
of UCMHs at a range of masses. While the previous strongest 
limits rely on the assumption that dark matter can annihi¬ 
late, the limits we find here are placed by gravitational meth¬ 
ods only, and are therefore equally applicable to any dark 
matter model. Our best limit of /mw 1.1 x 10“® is more 
than an order of magnitude better than previous ones. For 
masses 10“^ < M^/Mq < 10^^, our limits are the strongest 
available, even compared to earlier model-dependent ones. 

While different lensing methodologies appear to have 
been exhaustively used for investigating the small-scale 
structure of dark matter, time-delay methods have been 
mostly overlooked. This new methodology shows that time 
delays may yet provide the most sensitive measures of lens¬ 
ing to date. We have shown that, although small in ampli¬ 
tude, the gravitational time-delay signal due to UCMHs may 
indeed yet be seen in pulsar timing observations. Should a 
detection be made, these may be used to constrain mod- 
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els of early universe cosmology such as inflation and cosmic 
strings. 
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